Spectral statistics in an open parametric billiard system 
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We present experimental results on the eigenfrequency statistics of a superconducting, chaotic 
microwave billiard containing a rotatable obstacle. Deviations of the spectral fluctuations from pre- 
dictions based on Gaussian orthogonal ensembles of random matrices are found. They are explained 
by treating the billiard as an open scattering system in which microwave power is coupled in and 
out via antennas. To study the interaction of the quantum (or wave) system with its environment 
a highly sensitive parametric correlator is used. 
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Classical chaos manifests itself in universal spectral 
quantum fluctuations that can be described by random 
matrix theory (RMT) 0. While the earliest investiga- 
tions of spectral correlations were confined to nuclear 
physics 0, during the last twenty years the universality 
has been tested in other areas, like optical experiments 
, quantum dots Q , and acoustic setups || . The (local) 
spectral statistics depend generically only on the under- 
lying symmetries of the system. In particular, they are 
described by the Gaussian orthogonal ensemble (GOE) 
of real symmetric random matrices for spinless systems 
with time reversal symmetry, and by the Gaussian uni- 
tary ensemble (GUE) of complex Hermitian random ma- 
trices in the absence of time reversal invariance [|| . The 
sensitivity of the quantum (or wave) statistical properties 
to fundamental symmetries is obviously of great inter- 
est. For instance, it has been utilized to derive an upper 
bound for the magnitude of the time or parity violating 
component in nuclear interactions 

We investigate here spectral properties of a supercon- 
ducting microwave resonator where currents are induced 
by the measurement process. Although we study a specific 
wave system, the results are expected to be of general va- 
lidity in the physics of complex quantum systems (atoms, 
molecules, nuclei, quantum dots, ...). The influence of the 
flux of microwave power flowing from the feeding to the 
receiving antenna on the spectral properties of the system 
is so weak that it can only be detected through a highly 
sensitive diagnosis tool, in our case a parametric statisti- 
cal measure. In a previous experiment, the wave system 
was realized by a normal conducting microwave resonator 
attached to a large number of antennas [9j- There, the 
distribution of wave functions showed significant devia- 
tions from the GOE predictions, which were attributed 
to the transformation of the standing waves inside the 
closed microwave billiard into waves propagating from 
an emitting antenna into a large number of exit channels 
|lflj . The aim of the present paper is to go further into 
the investigation of this mechanism. We will show that 
deviations from GOE behaviour are already observed in 
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FIG. 1: Sketch of the billiard, showing the outer boundary, as 
well as the rotatable wedge-shaped piece of teflon. The angle 
of rotation a of the wedge defines the parameter (its initial 
orientation is arbitrary). 



a resonator with only three (or less) attached antennas, 
when studying spectral properties as a function of a pa- 
rameter. 

The experiment discussed here has been performed 
with a superconducting microwave resonator, whose 
high-quality factor is typically Q — 10 5 or larger [i~l|| . 
i.e. dissipative processes in the resonator are reduced 
to a minimum, thereby ensuring a high spectral resolu- 
tion. Results obtained with a flat cylindric resonator are 
presented. Aside their intrinsic interest, such resonators 
mimic two-dimensional quantum billiards of correspond- 
ing shapes [III UJI LU1 LLii ■ The analogy is based on the 
isomorphism between the scalar Helmholtz equation of 
the electric field E for wavelengths longer than twice the 
height of the resonator, where E(r) = ^(x,y)e z is per- 
pendicular to the billiard (xy) plane, and the Schrodinger 
equation for the wavefunctions in the quantum billiard. 
The eigenvalues kf of a closed resonator satisfy the 
Helmholtz equation with Dirichlet boundary conditions 
imposed on *S>(x,y). They are directly related to the 
eigenenergies of the corresponding quantum billiard. In 
this analogy, the Poynting vector plavs the role of the 
quantum probability current density 15]. 

The microwave resonator has been manufactured from 
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FIG. 2: Part (# 200- # 250) of the experimental eigenvalues 
on an unfolded, i.e. dimensionless scale as a function of the 
parameter a. 

lead-plated copper, as in ^(|; its shape is shown in Fig.^ 
During the measurements it has been placed in a liq- 
uid helium cryostat at a temperature of T = 4.2 K, 
which guarantees superconductivity of the lead surface. 
The outer boundary of the resonator has the shape of 
a desymmetrized straight-cut circle. The dynamics in- 
side the corresponding classical billiard is chaotic 0] . A 
dielectric wedge of teflon inside the resonator has been 
rotated with a leverage from outside the cryostat. The 
cavity has been coupled to one feeding and two receiving 
antennas. They are tiny metal pins of 0.5 mm in di- 
ameter and have been adjusted such that they mechani- 
cally reach only some 100 microns into the interior of the 
cavity guaranteeing weak coupling and hence minimal 
disturbance of the excited field in the resonator. Using 
a HP-8510C network analyzer the transmission spectra 
have been measured in the frequency range up to 18 GHz 
for 37 equidistant settings of the angle a (cf Fig. Q in 
steps of 2.5 degrees. 

Complete eigenvalue sequences of 440 resonances for 
each value of the parameter a have been measured. Fig- 
urc[2]shows the eigenvalues number 200 to 250 as a func- 
tion of the parameter, the so-called eigenvalue dynam- 
ics or parametric fluctuations, where for each of the 37 
spectra the eigenvalues have been unfolded by scaling 
them to unit mean spacing [ill Il2^ . As the parame- 
ter a is varied, the positions of the resonance frequen- 
cies (i.e. the real part of the resonance) describe irregu- 
lar oscillatory curves that generically do not cross each 
other. The oscillations have a mechanical interpretation, 
namely the derivative —dkf/da is proportional to the 
torque exerted by the i-th electromagnetic eigenmode on 
the teflon wedge. For the statistical analysis of the spec- 
tra the whole set of resonance frequencies has also been 
unfolded with respect with respect to the parameter a by 
following the proceedure described in . This allows to 
properly incorporate the characteristic scales associated 
with the frequency and parameter secular variations 0] , 
thereby defining dimensionless quantities. 

In Fig. [31 the spacing distribution P(s) of the distance 




FIG. 3: The experimental nearest-neighbor spacing distribu- 
tion (histogram); the GOE random matrix result (full line). 
The eigenvalues # 331 to # 440 of each of the 37 spectra have 
been used. 

between consecutive eigenvalues computed at fixed val- 
ues of a is shown. The experimental P(s) is in good 
agreement with the GOE result. A similar agreement 
with GOE is found for the number variance E 2 and the 
least mean square statistics A3. This is in agreement 
with common expectations, especially in the present ex- 
periment with a high Q value of the resonator and a 
small coupling to the antennas. However, a minimum 
coupling is unavoidable (open system), thereby implying 
a presumably small perturbation of the closed system. 
The purpose of the present investigation is the identifi- 
cation of signatures of this disturbance produced by the 
measurement process in the spectral properties. 

It is by now well established that in systems depend- 
ing on parameters, correlations between eigenvalues at 
different parameter values lead to important extensions 
of RMT universalities 0, |2jj. In the present exper- 
iment, these correspond to correlations between spec- 
tra at different orientations of the teflon wedge. Several 
parametric correlators \2l\ have been computed, namely 
the velocity distribution, the curvature distribution, the 
velocity-velocity correlator, the diffusion correlator and 
the distribution of the spacings at avoided crossings (a 
general presentation of the results will be given else- 
where). We focus here on the latter, i.e. the probability 
distribution P(c) of the local minima c of the distance 
between neighboring levels as the parameter a is varied. 

The parameter dependence of the chaotic resonator is 
modeled by the following ensemble, 

H(p) = cos fi ■ Hq + sin fi ■ Hi , (1) 

where Hq, H\ are iV-dimensional GOE random matrices 
and fj, is a real parameter. Before performing statistical 
analyses, the spectra are unfolded (with respect to energy 
and parameter dependence, see above). We will denote 
by [i resc the resulting rescaled parameter. Though the 
form for large N of the probability distribution P(c) has 
not been derived for this model, it is well approximated 
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FIG. 4: The experimental distribution of avoided crossings 
(full line histogram) compared to the GOE prediction Eq.© 
(full line curve). The dashed line histogram is obtained from 
the model (JIJ with A = 0.02 and using a discrete set of 
rescaled parameters firesc (see main text). The eigenvalues 
# 331 to # 440 of each of the 37 spectra have been used. 

by the N = 2 result jH 

2 c 2 
P GOE (c)«-exp( ), (2) 

where the scale of c has been chosen such that (c) = 1. 

We have evaluated P(c) from the experimental data. 
The result is presented in Fig. It is compared to the 
prediction Eq. Though the general trend shows an 
overall agreement, systematic deviations are observed. In 
particular, a lack of small spacings of avoided crossings 
is clearly visible in the experimental distribution. 

It is interesting to note that similar results were ob- 
tained in |23| . where they were attributed to an insuffi- 
cient experimental resolution in the parameter variation. 
This was also our first suspicion and we tested it nu- 
merically with a random matrix model. The rotatable 
obstacle, i.e. the varying of the shape of the resonator, is 
accounted for in the model by considering the parameter 
dependent ensemble of random matrices given in Eq. 
where the dimension of the matrices Hq, H\ is chosen as 
iV = 1000. Only the 300 central eigenvalues of each di- 
agonalization are used. If the rescaled parameter fj, resc is 
treated as a continuous parameter then a curve for P(c) 
very similar to Eq.J2J) results. If, in contrast, [i resc is dis- 
cretized such that a set of parameter values very similar 
to the experimental one is used, then the agreement with 
the experimental histogram improves. In particular, a 
dip at c ~ appears. However, systematic deviations of 
the experimental P{c) and also of other parametric statis- 
tical measures studied (velocity and curvature distribu- 
tion, velocity-velocity correlator) from the corresponding 
predictions of the model (JIJ persist. We thus conclude 
that the disagreement between the RMT predictions of 
the model Q and the experimental results cannot be 
entirely attributed to an insufficient resolution with re- 
spect to the parameter. Other effects, like the stability 
of the dip with increasing frequency, were also checked. 



We have devided the 400 levels in windows of increasing 
frequency and computed P(c) for each of them. We were 
not able to detect any systematic trend in the behavior 
of the dip as frequency increases. The present available 
data thus exclude the attribution of the dip to a finite 
size effect. 

Since the experiment is performed with a supercon- 
ducting resonator, in the absence of the antennas and 
neglecting dissipation, the cavity is an isolated time- 
reversal invariant system which, in principle, should be 
correctly described by the parametric statistical model 
Eq. (JTJ. However, for the measurement of a spectrum, a 
typical procedure is to couple the system to the exterior 
through antennas, and to emit an input signal via one 
of them and receive the output signal via another one. 
Hence, the effective Hamiltonian describing the spectral 
properties of the resonator is non-Hermitian (open sys- 
tem), and the statistical properties of the spectrum are 
not expected to coincide with those predicted from the 
Hermitian model Eq. . If the coupling to the antennas 
is weak, a small but nonzero change in the position of the 
real part of the resonances with respect to the closed sys- 
tem is thus expected. As Fig. shows, with the present 
experimental conditions this shift in the resonance fre- 
quencies has no visible effect on the spectral fluctuations 
24] for a fixed value of the parameter. In contrast, as 
we will see below, it induces sizeable deviations in the 
parametric statistical properties. 

In the present experiment the antennas act as single 
scattering channels as their diameter is small compared 
to the wavelengths of the microwaves in the total fre- 
quency range. Wave scattering in such a three-port sys- 
tem is described by a 3 x 3 scattering matrix of the form 

m mum 

S = I + 2iW T (H(fi)-iWW T - Eiy'tw , (3) 

whose derivation is based on the theory of quantum 
scattering (formulated e.g. in [29j). Here, / is the 
identity matrix and H(/j,) the Hamiltonian of the res- 
onator. It is modeled by the N x N parametric GOE 
defined in Eq. The matrix W is an N x 3 matrix, 
W oc (Xi,X2,X 3 ), that describes the coupling of the 
resonator to the antennas (the j-th component of the N- 
dimensional column vector X a couples the j-th internal 
wavefunction to the tr-th antenna) . From Eq. @ , the res- 
onances are obtained as the eig envalues of the effective 
non-Hermitian Hamiltonian |3ij 

H eff (X, tj) = H{fj)-i\^= (XxX? + X 2 X% + 1 3 A 3 T ) . 

(4) 

Since the system is time reversal symmetric and the cou- 
pling is weak, the emission of waves from one antenna 
and its detection in another is modeled with real col- 
umn vectors X a as in |2J| . Consistently with the random 
model adopted for the Hamiltonian of the resonator, they 
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are considered as independent random variables with a 
Gaussian distribution whose width is set to unity. Then, 
in Eq. Q the parameter A measures the strength of 
the coupling of the resonator to the antennas in units of 
the mean spacing tt/VN of the eigenvalues of £f(/i). In 
contrast to the present experiment, the strong coupling 
regime has been investigated in [2(| |^| . 

We have studied numerically the statistical properties 
of the eigenvalues of H e f f . For small values of A the res- 
onances are close to the real axis and tend to the eigen- 
values of H(fi). As A increases, and up to A ~ 0.5, the 
imaginary part of the resonances increases. For larger 
values, the resonances split into two groups: three of the 
resonances, their number corresponding to the rank of 
the perturbation, move deeply into the complex plane 
while the remaining N — 3 approach again the real axis 
with increasing A. Due to the weak coupling, the present 
experiment should correspond to relatively small values 
of A. We find numerically that for values of A smaller than 
« 0.05 the ratio of the imaginary to the real part of the 
eigenvalues of H e ff is smaller than 0.005, and its ratio to 
the mean spacing between adjacent real parts is less than 
0.1. The numerical P(s), the A3-statistics and the S 2 - 
statistics agree with GOE, in accordance with the exper- 
imental results. The distribution P(c) of avoided cross- 
ings, however, deviates from the model for A larger 
than about 0.01. Before reaching again a GOE-like be- 
havior (J2J at A ~ 2.5 as tested numerically, a sharp <5-like 
peak at c ~ followed by a dip is observed. The peak 
size increases until A ~ 0.5, and then decreases. This be- 
havior differs both from the GOE behavior, Eq.(J2|, and 
from the experimentally observed distributions. A more 
careful analysis of the behaviour of the real and the imag- 
inary part of the eigenvalues of H e ff, Eq. Q, as a func- 
tion of the parameter fi shows that the contributions to 
the peak of P(c) at c = are due to crossings of the real 
parts of two complex eigenvalues. Incidentally, this be- 
haviour is characteristic of non-Hermitian Hamiltonians 
studied recently j^- 

Why is the sharp peak at c ~ predicted by the model 
(|3"|) not observed in the experiment? Our interpretation 
is that this is due to the discrete variation of the ex- 
perimental parameter. Because of the discrete sampling, 
the probability to observe small spacings is strongly re- 
duced. And indeed, the smallest spacing observed in the 
unfolded experimental data is about s ~ 0.014. To com- 
pare theory and experiment a discretization of the pa- 
rameter in the model (@J has to be performed. In Fig. 01 
such a comparison is made using in Eq.(@J a discrete 
set of parameter values \x resc whose step size is similar 
to the experimental one, and A = 0.02. A good over- 
all agreement between both curves is obtained, not only 
around the dip close to the origin. A similar agreement 
is obtained for the other parametric correlators studied 
(curvature distribution, velocity- velocity correlator, etc) , 
thus providing a globally consistent picture of the exper- 



imental data. 

The experimental as well as the theoretical results thus 
indicate that, while absent in the nearest neighbor spac- 
ing distribution, we are observing in more sensitive para- 
metric spectral functions signatures induced by the mea- 
surement process. Interestingly, similar mechanisms were 
recently studied in the physics of cold atoms 33]. Two 
ingredients are important in order to understand the ex- 
perimental results. First, we had to model the system 
by incorporating its coupling to the external world (flux 
is fed into the resonator and coupled out via one or two 
antennas). Second, we had to take into account the dis- 
creteness of the parameter variations. 

The proposed parametric " spectral detector" is clearly 
a powerful tool to study the interaction of a quantum 
(or wave) system with its environment. The underly- 
ing working principle is very general since it depends on 
a fundamental physical principle, namely the standing 
waves inside the closed resonator are transformed into 
waves propagating from an entrance antenna to an exit 
antenna. In particular, it may be useful in the analysis of 
more controlled experiments concerning the interplay be- 
tween the measurement process, the currents it induces 
through the cavity, and the dissipative processes. 
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